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Abstract. We present new imaging and spectroscopy 
of the post-red supergiant binary AFGL 4106. Corono- 
graphic imaging in Ha reveals the shape and extent of the 
ionized region in the circumstellar envelope (CSE). Echelle 
spectroscopy with the slit covering almost the entire ex- 
tent of the CSE is used to derive the physical conditions 
in the ionized region and the optical depth of the dust 
contained within the CSE. 

The dust shell around AFGL 4106 is clumpy and mixed 
with ionized gas. Ha and [N ii] emission is brightest from 
a thin bow-shaped layer just outside of the detached dust 
shell. On-going mass loss is traced by [Ca ii] emission and 
blue-shifted absorption in lines of low-ionization species. 
A simple model is used to interpret the spatial distribution 
of the circumstellar extinction and the dust emission in a 
consistent way. 

Key words: binaries: spectroscopic — circumstellar mat- 
ter — Stars: individual: AFGL 4106 — Stars: mass loss 
— Stars: AGB and post-AGB — supergiants 



1. Introduction 

In their final stages of evolution, both massive (Af^8 Mq) 
and intermediate-mass (iSSAfSSS M0) stars can become 
dust-enshrouded as a result of intense mass loss at rates 
of 10^^ up to 10^^ M0 yr^^. For the massive stars, this 
heavy mass loss occurs during the red supergiant (RSG) 
phase, while for the intermediate-mass stars this happens 
during the Asymptotic Giant Branch (AGB) evolution. 
After having lost a significant fraction of their initial mass 
during this episode, the photospheric temperature of the 
star increases again and the stellar ultraviolet (UV) radia- 
tion field starts ionizing the dusty CSE, producing either a 
Planetary Nebula (FN) or a post-RSG nebula surround- 
ing a Wolf-Rayet star. The mechanisms at play during 
the short-lasting transition stage between RSG/ AGB and 
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WR/PN are still poorly known. Especially the final evolu- 
tion of the massive objects is hitherto poorly documented. 
One such massive object is AFGL 4106. 

Garcfa-Lario et al. (1994) discovered nebular line emis- 
sion from the IR object AFGL 4106, which came as a sur- 
prise because the central star was thought to be of spectral 
type G: too cool to ionize a CSE. The dust shell around 
AFGL 4106 was imaged in the infrared (IR) by Molster 
et al. (1999) who modelled the spectral energy distribu- 
tion as observed by ISO with a radiation transfer code. 
They provide evidence that the object is a binary consist- 
ing of a late-A/early-F type star and a somewhat fainter 
M-type companion. Their distance estimate yields lumi- 
nosities too high for AGB evolution and hence AFGL 4106 
must be the result of post-RSG evolution. 

Here we present and analyse new observations of this 
object and focus on the spatial distribution and physical 
conditions of the CSE around AFGL4106. Coronographic 
imaging in Ha shows the spatial distribution of the ionized 
material around AFGL 4106. Spatially resolved echelle 
spectroscopy is performed to measure expansion veloci- 
ties, electron densities and the internal extinction by the 
dust. The spatial distribution of dust extinction and emis- 
sion is modelled in a consistent way. Finally we discuss 
the structure of the circumstellar envelope and the recent 
mass-loss history of AFGL 4106. 

2. Observations 

2.1. Ha coronographic imaging 

On February 4, 1996, we used the multi-mode instrument 
EMMI at the ESO 3.5m NTT on La Silla, Chile, to image 
AFGL 4106 through a filter centred at Ha, with a rather 
narrow width of 33 A to avoid contamination by [N li] 
emission. Because of the high apparent brightness of the 
star (9*'' mag in V) coronographic techniques were used to 
limit the saturation of the CCD. A glass plate was inserted 
in the aperture wheel, holding six dark blots of different 
size and attenuation. The star was placed behind a 1.9" 
blot with a central attenuation of a few mag. Two frames 
of each 5 min integration time in Ha were combined. A 
10 s exposure through a broad-band Johnson R filter was 
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used to correct for the continuum contribution within the 
Ha filter band. The R-band image was aligned with the 
Ha image and scaled linearly by comparing the intensi- 
ties of several field stars. Before subtracting the R-band 
image from the Ha image both were corrected for struc- 
ture in the instrumental response over the pixels of the 
CCD by dividing them by images of the morning twilight 
sky, taken through the corresponding filters. Remaining 
CCD artifacts and cosmic ray impacts were removed by 
hand. The seeing was constant at 1", and the pixel size 
was 0.268". 

2.2. Echelle spectroscopic imaging 

On February 4, 1996, we used the same telescope and 
instrument to take cchcUc spectra of AFGL 4106 with 
several slit orientations. The wavelength region extended 
from 6000 to 8350 A. Using a slit width of 1" a resolv- 
ing power of i? ^ 7.5 x 10** was obtained corresponding 
to a velocity resolution of 4 km s~^. The slit length was 
set to 15" and the spatial resolution was 0.268" per pixel. 
We centred the star in the middle of the slit but obtained 
15 min spectra with different position angles: 0°, 45°, 90° 
and 135° where the position angle is defined from West 
over North. The reduction of the spectra was standard 
and included flat-fielding, wavelength calibration on the 
basis of Th-Ar exposures, spectral response correction on 
the basis of a measured spectrum of standard star HD 
60753, interactive cosmic hit cleaning, extinction correc- 
tion with average extinction values, and finally absolute 
flux calibration. The seeing was about 0.8 to 0.9". 

3. Results 

3.1. Ha coronographic imaging 

The coronographic image in Ha is displayed on a loga- 
rithmic flux density scale in Fig. 1. North is up, and East 
is to the left. The star itself is heavily saturated despite 
the use of the coronographic blot, resulting in the artificial 
NS-orientated bar. The absolute fiux density calibration 
was derived from the spectra and is estimated to be accu- 
rate within ^ 20%. The flux density level of the brightest 
part of the extended emission is ~ 7 x 10~^^ W m~^/zm~^ 
arcsec"^. The faintest emission visible in the picture is 
about a factor of ten fainter. 

AFGL 4106 appears to be situated in a bright, spa- 
tially extended emission complex. From N to W of the star, 
the strongest emission delineates a bow-shaped structure 
much akin the bow-shocks associated with stars that move 
supersonically through the interstellar medium (Kaper et 
al. 1997). The emission bow is located at a projected dis- 
tance of 5" up to 10" from the central star. The faintest 
Ha emission is detected up to ^ 14" from the central 
star. In the SE much less Ha emission is detected. Inter- 
estingly, the 10 /xm emission (Molster et al. 1999) shows 
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Fig. 1. Coronographic image in the light of Ha on a log- 
arithmic flux density scale. The flux density level of the 
brightest part of the extended emission is~7xlO^^^W 
m~^/zm~^ arcsec^^. The faintest emission visible in the 
picture is a factor of ten fainter. North is up, and East is 
to the left. 

a clear anti-correlation with the Ha emission: the 10 fim 
emission peaks in the SE and is faintest in the NW. 

3.2. Echelle spectroscopic imaging 
3.2.1. Absorption lines 

The spectral type has been redetermined by van Winckel 
et al. (in preparation) and Molster et al. (1999). The opti- 
cal spectrum is dominated by a late-A or early-F type star 
(Teff ^ 7500 K) with a nitrogen-enhanced photosphere 
resulting in strong N I absorption lines around 8000 A. 
They also spectroscopically discovered the presence of an 
early-M type companion star (Teff ^ 3750 K). The self- 
absorbed Ha line profile centred on the star is identical to 
that observed by Garcfa-Lario et al. (1994). This suggests 
that the Ha emission conditions in the vicinity of the star 
have been stable over a period of at least six years. 

A total of 31 absorption lines of the atomic and singly 
ionized species Fe i, Fe ii, Sc ii, Si ii, O i, Mg ii, Na i and 
N I were selected for measuring the star's radial velocity. 
Considering a heliocentric correction for the movement of 
the Earth of -1-13.3 km s~^, the heliocentric velocity of the 
star is determined at v^, = +8.1 ±0.7 km s"'^. For an indi- 
vidual line, the typical deviation from the mean was 4 km 
s~^, which is equal to the spectral resolution. The stellar 
velocity differs considerably from the value of ?,v = — 47±2 
km s^^ derived by Garcia-Lario et al. (1994) from the Si ii 
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Fig. 2. Radial velocity curve for AFGL 4106. The solution 

is not unique, but periods mueli longer than 12 yr are 
difficult to reconcile with this data. 
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Fig. 3. Spectrum around Si ii A6347 A (sohd) and the 
same spectrum mirrored with respect to 6347.29 A where 
the line profile is deepest (dashed). Below, their difference 
spectrum is plotted. The wavelength axis is heliocentric. 



AA6347,6371 lines (after transforming from LSR to hclio- 
centre) in their spectra taken in February 1990. In March 
1992, ^^CO J=1^0 emission at 115 GHz was detected 
by Josselin et al. (1998). From the mean of the veloci- 
ties of the blue- and red-most CO emission we estimate 
Vi, = — 3±1 km s~^. An additional measurement from pho- 
tosphcric absorption lines in high resolution spectra taken 
in July 1996 (van Winckel et al., in preparation) yields 
Vi, = —4 ± 1.2 km s~^. An approximate radial velocity 
curve (Fig. 2) for the maximum period possible, assuming 
circular orbits, yields Pmax = 4500 d. For (currently) equal 
masses of the two stars Kepler's third law yields the mass 
of each star: M ^ 15sin^'^i M0 with orbital inclination 
angle i (shorter periods imply smaller masses). 

The Si II AA6347,6371, Fe 11 A6456, N i A7468 and Fe i 
A8327 absorption lines have asymmetric line profiles, with 
more absorption in the blue wing. This may be caused 
by outflowing matter due to present day mass loss. The 
spectrum around Si 11 A6347 is shown as an example (Fig. 
3). The spectrum was normalised by dividing by a con- 
stant equal to the spectrum level outside of the absorption 
line, i.e. the spectral slope is not affected. Also plotted is 
the difference between this spectrum, and the spectrum 
after mirroring with respect to the absorption extremum 
at 6347.29 A(after heliocentric correction). This difference 
spectrum shows the absorption excess in the blue wing of 
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Fig. 4. Spatial and kinematic map of [N 11] A6583. The 
stellar continuum is subtracted, and the flux scale is loga- 
rithmic. The relative positions are defined as running from 
E to W at a position angle 0°, SE to NW at 45°, S to N 
at 90° and SW to NE at 135°. The wavelength axes are 
transformed into heliocentric Doppler velocity axes. 



the line profile, reaching a maximum at a blueshift of ^ 40 
km s~^. 



3.2.2. Diffuse IntersteUar Bands and the E{B - V) 

The centroid wavelengths and equivalent widths of 68 
Diffuse Interstellar Bands (DIBs) were measured. These 
DIBs are found over the entire spectral coverage, although 
few are found at wavelengths longer than 7400 A. The 
DIBs were identified using Jenniskens & Desert (1994), 
from which also the conversion factors between equiva- 
lent width and colour excess E{B — V) for each DIB were 
adopted. The derived colour excess is E{B — V)= 0.9t.Q l 
mag, confirming the estimated E{B — V) = 1.0 ± 0.2 mag 
from a smaller selection of DIBs by Molster et al. (1999) to 
whom we refer for a thorough discussion on the distance 
to AFGL 4106. The heUocentric velocity is vbib = +3±2 
km s^^. The velocity of the DIBs differs significantly from 
both the velocity of AFGL 4106 as measured from the 
spectra as well as the system velocity of AFGL4106, con- 
firming that the DIBs are of interstellar not circumstellar 
origin. 
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3.2.3. [N ll] emission lines 

The [N ii] A6583 line is the strongest emission line in 

AFGL 4106, and the best tracer of the kinematics of the 
circumstellar nebula. Unlike the Ha, there is no stellar 
contribution to the line profile. The underlying stellar con- 
tinuum is subtracted by scaling the spatial profile of the 
stellar spectrum outside the kinematic extend of the [N ii] 
emission. In this way a map is constructed of the [N ll] 
emission as a function of the position along the slit relative 
to the stellar position, and the heliocentric Doppler veloc- 
ity, for each of the four slit orientations (Fig. 4). The flux 
scale is logarithmic, and the relative positions are counted 
from E to W (at 0°), from SE to NW (at 45°), from S to 
N (at 90°) and from SW to NE (at 135°). Close to the 
stellar position the stellar continuum was too dominant 
over the [N ii] emission to reliably correct for it, and this 
region in the four maps is covered by the labels with the 
slit position angles. 

The maps of the [N ii] emission are consistent with 
emission from a radially expanding shell, with an expan- 
sion velocity of ^ 30 km s^^ and a radius of ~ 7". As in 
the Ha emission map, the NW side is brighter than the 
SE side. The expansion velocity is smaller than the 40 km 
s~"^ derived from the Si ii absorption and also smaller than 
the Wexp = 35 ± 1 km s~^ that we estimate from the ^^CO 
J=l— >0 emission at 115 GHz detected by Josselin et al. 
(1998). 

The spectra at the four different position angles were 
summed and integrated over the blue- and red-shifted ve- 
locity range, respectively. This results in an effective cross 
section of the spatially extended emission complex run- 
ning roughly from SE to NW, i.e. from the minimum to 
the maximum of the emission measure. The result is plot- 
ted in the top panel of Fig. 5. Assuming that the source of 
the [N ll] emission is a spherically symmetric expanding 
shell, the logarithm of the ratio of the blue- over the red- 
shifted emission corresponds to the optical depth within 
the emitting shell. This is because the red-shifted emission 
corresponds to the far side of the shell and has to traverse 
through the dusty CSE within the shell before it emerges 
at the observer's side which is also where the blue-shifted 
emission arises from. The derived spatial profile of the op- 
tical depth is shown in the bottom panel of Fig. 5. 

The optical depth profile traces the dusty CSE out to 
a radius of 5", which corresponds to the spatial extent 
of the 10 /um emission (Molster et al. 1999), to the inner 
boimdary of the Ha maximum (Fig. 1) and the forbidden 
line maxima (Fig. 4). The optical depth at A = 6583 A 
exceeds unity at a distance of ^ 2" from the star, but is 
lower closer to the star. This corresponds to the inner cav- 
ity also seen in the 10 /xm emission. The red-shifted [N ii] 
emission becomes stronger than the blue-shifted emission 
at relative positions ^5". This cannot be explained by op- 
tical depth effects and must therefore indicate deviations 
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Fig. 5. Spatial profile of [N ii] A6583 (top panel), obtained 

by summing the spectra taken at the four different posi- 
tion angles, with the relative position defined as in Fig. 4. 
Positive relative positions roughly correspond to NW. The 
blue- and red-shifted emission is shown individually. The 
logarithm of their ratio yields the internal optical depth 
r (bottom panel) under the assumption that the [N ii] 
emission arises from a spherically symmetric expanding 
shell. 

from spherical symmetry in the density and/or excitation 
conditions in that part of the nebula. 

3.2.4. [S ii] emission lines 

The [S ii] AA6717,6731 emission is co-spatial with the [N ii] 
emission, as seen from the spectra at the four different 
position angles. The spectra arc summed in the same way 
as in preparing the [N ii] spatial profile in Fig. 5, but now 
the emission line is integrated over its entire kinematic 
extent. This results in one-dimensional spatial emission 
profiles for both components of the [S ii] doublet (Fig. 
6, top panel). The emission peaks in the NW (positive 
relative position), where the Ha is brightest too. 

The ratio of the emission profiles of the 6717 and 6731 
A components (Fig. 6, bottom panel) may be used to de- 
rive the local electron density rig (Ostcrbrock 1988). For 
Ue < 100 cm~^ this ratio is > 1.3, with a maximum of 1.42. 
For Tie > 10^ cm~^ the ratio is < 0.5, with a minimum of 
0.44. This holds for an electron temperature of = 10^ 
K. The derived electron density scales with y/T^. The CSE 
of AFGL 4106 outside of the region of 10 /zm emission has 
rig < 100 cm~"^, whereas the inner CSE corresponding to 
the cavity in the 10 /um emission has > 10** cm~'^. 
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Fig. 6. Spatial profile of [S ii] A6717 and A6731 (top 
panel), obtained as in Fig. 5 but integrated over the entire 
kinematic extent of the emission line. Their flux ratio is 
plotted in the bottom panel. 

3.2.5. Ha emission line 

Spatial profiles of the Ha emission were created in the 

same way as for [S ii]. The Ha profile is compared with 
the profiles of [N ii] (Fig. 7) and [S ii] (Fig. 8). 

The flux ratios of Ha and forbidden lines are 
in the range log(FHa/-F'[N ii]) ^ —0.3 to —0.05 and 
log(i^Ha/-F[s II]) ~ 0.3 to 0.7. This strongly suggests ex- 
citation conditions similar to those in Planetary Nebulae 
but different from those in H II regions (Garcia-Lario ct 
al. 1991). The Ha emission is co-spatial with the forbid- 
den line emission. The Ha/[N ii] ratio seems smallest in 
the brightest part of the nebula around a relative position 
of +5". The Ha/[S ii] ratio seems to be smallest close to 
the star. This effect may be caused by the spatial varia- 
tion of the ratio of the [S ii] components, with the A6717 
component essentially following the Ha emission. 

3.2.6. Fe I and [Ca ii] emission lines 

Though reported by Garcia-Lario et al. (1994), 16.6 A 
wide Fe I emission centred at 6380.7 A could not be con- 
firmed. 

Instead, weak fluorescence emission of the IF doublet 

of [Ca ii] AA7291,7324 is detected (Fig. 9). The spectra 
have been normalised in the same way as the Si ii spec- 
trum (Fig. 3). The emission is centred on the star and 
spatially unresolved. The line profiles are distorted due 
to the presence of numerous telluric absorption lines but 
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Fig. 7. Spatial profiles of Ha and [N ii] AA6548.6583, ob- 
tained as in Fig. 5 (top panel). Their fiux ratio is plotted 
in the bottom panel. 
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Fig. 8. Spatial profiles of Ha and [S ii] AA6717,6731, ob- 
tained as in Fig. 5 (top panel). Their flux ratio is plotted 
in the bottom panel. 



appear to be kinematically symmetric with respect to the 
stellar velocity. The emission extends between about ±20 

km s""'^, which is somewhat narrower than the Ha, [N ii] 
and [S ii] emission. The equivalent widths of the lines are 
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Fig. 9. Spectrum around [Ca ii] AA7291,7324. The wave- 
length axes are transformed into heliocentric Doppler ve- 
locity axes. 

M^729i = —26 ± 3 and VF7324 = —31 ± 4 mA, respectively. 
The line profiles peak at a flux density (above the stellar 
continuum) of F7291 ~ 3.0 x IQ-^^ ^nd ^^7334 ~ 2.3 x IQ-^^ 
W m~^/xm~^, respectively. 

The location of the [Ca 11] emission can be further con- 
strained by estimating where the outflow reaches the crit- 
ical density for the [Ca II] lines of ricrit = lO'^ cm~^. How- 
ever, this requires c;alculation of the mass-loss rate. The 
mass-loss rate may be estimated from the CO emission 
reported by Josselin et al. (1998). For the moment it is as- 
sumed that this mass-loss rate is valid for the unresolved 
inner regions of the nebula. 

With distance Z?, outflow velocity v, CO abundance / 
with respect to H2, CO peak antenna temperature T, and 
a simplified model for the CSE (Kastner 1992), the total 
mass-loss rate can be derived: 



-12 



Kom) 



M = 3.8 X 10^ 

X ^Moyr-i 

with the parameter 7 approximated as follows: 



\ for M < 10-^ Mq yr 
7 = I for 10""^ <M < 10-6 IVI0 yr 
f for 10-6 <M < 10-5 M, 



(1) 



(2) 



For higher mass-loss rates, more detailed modelling is re- 
quired. In that case the CO line profile will become sat- 
urated, and hence it is no longer the peak antenna tem- 
perature, but the shape of the line profile that yields an 



estimate of the mass-loss rate. With v - 
/ = 2 X lO-'', T = 0.14 K, and 7 = f , 
rate is obtained: 

M = 3.5 X 10-^£>^ [kpc]^ Mgyr-^ 



= 35 km s ^, 
the mass-loss 



(3) 



This implies a very high mass-loss rate if the star were 
to be considerably more distant than one kpc, but still an 
order of magnitude less than the mass-loss rate derived by 
IVEolster et al. (1999) who estimate M ~ 3 x IQ-^D [kpc] 
M0 yr-i. 

The angular separation d from the star from which the 
[Ca 11] emission arises can now be estimated assuming a 
constant outflow velocity of w = 35 km s-^, and an average 
particle mass of 2.7 x 10 



-24 



d = 1.4 X lO^^D [kpc] m = 0.09" 



(4) 



Hence for distances of D^20 kpc the [Ca 11] emission must 
arise from well within the inner cavity that is seen in the 
10 /xm emission and in the internal extinction of the [N ll] 
emission. If the mass-loss rate in those inner regions is 
lower than the CO estimate, then the [Ca 11] emission is 
located even closer to the star. This is consistent with the 
[Ca 11] emission being spatially unresolved. 

4. Modelling the circumstellar dust envelope 

The dust around AFGL 4106 causes extinction of line 
emission that originates from behind the dusty CSE, and 
emission in the mid- and far-IR. Our spatially resolved 
spectra of the nebular line emission were used to construct 
a spatial profile of the optical depth of the dusty CSE, and 
the spatially resolved dust emission at 10 /xm was pub- 
lished by IVEolster et al. (1999). Here we make an attempt 
to simultaneously model the spatial profiles of the optical 
depth and of the dust emission. 

4.1. Circumstellar dust extinction 

Suppose that the bulk of the Ha and [N 11] emission orig- 
inates from a region just outside the dusty envelope. The 
emission from the front and backside is observed at blue- 
and red-shifted velocities, respectively. Assuming equally 
strong front and backside emission, the ratio of the ob- 
served line flux at red- and at blue-shifted velocities di- 
rectly yields the differential extinction experienced by the 
backside emission after traversing the dusty CSE. The 
[N 11] A6583 line is best suited for this method because 
it is strong and its intrinsic width is small compared to 
the much larger thermal width of the Ha line emission. 

The mass-loss rate history through the dusty CSE is 
parameterised as: 



Mr 



(5) 



The case of a = corresponds to a steady mass loss. The 
outflow velocity through the dusty CSE is assumed to be 
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Fig. 10. Schematic representation of the dust shell around 
AFGL 4106, and the coordinates used in the modelling. 

constant. The optical depth ro.ee at 0.66 /xm as a function 
of dimensionless projected distance A is: 



TO.66 = Cfi'(A) 



(6) 



where A = 1 at the outer radius R of the dusty CSE, and 
A = Ao at the inner radius Rq of the dusty CSE (Fig. 
10). The spatial profile of the observed optical depth has 
a shape: 



5(A) 



J X 



Fig. 11. Top panel: spatial profile of the optical depth 
r at 0.66 jim as derived from the [N il] A6583 emission. 
Bottom panel: spatial profile of the 10 jim emission, for 
p — 1. In both panels, calculated profiles are plotted for 
different mass-loss rate histories (a, see text). 

4-.2. Circumstellar dust emission 

The spatial profile of the observed N-band flux density per 
square arcsecond can be calculated as follows: 



with lower integration boundary: 

1 if A > 1 
A if Ao < A < 1 
Ao if A < Ao 

and a scaling parameter: 



c 



^i^MR ^ Qo.66 



(7) 



(8) 



(9) 



Fiv = C'5'(A) 
with a shape: 

5'(A) = 



y^ 



f{y)dy 



_ ^2 

and a scaling factor: 

1 V 2hc^ Qn 



SttRvs a 

with dust-to-gas ratio tp = 5 x 10~^, constant outflow 
velocity v, dust grain specific mass s = 3 gcm~^ and 
radius a, and extinction coefficient Qo.ee = 14(i [/um] for 
oxygen-rich dust (Volk & Kwok 1988). 

For a e {—2,-1,0,1,2} the function g can be eval- 
uated analytically. To obtain a reasonable fit to the ob- 
served profile, the parameters Aq, R and C. need to be m = [e 
tuned accordingly. The model profiles are convolved with 
a gaussian of ct = 0.36", corresponding to the Full Width ^ 
at Half Maximum of the stellar continuum in the spatial 
direction (0.85"). We fitted by eye, aiming at equal max- 
imum optical depth for all a. The results are plotted in 
Fig. 11 (top panel) and summarised in Table 1. 



648, 000 



^% Qo.66 



c 



(10) 



(11) 



(12) 



The wavelength dependence of the extinction coefficient Q 
for oxygen-rich dust (Volk & Kwok 1988) is convolved with 
the TIMMI N-band filter curve, yielding Qat = 3.18a [^m]. 
The dust emission is given by a blackbody and experiences 
extinction on its way out of the CSE: 



^^P ( A,vfcT(j/) ) 



g-TN+Tiv(l/) 
^ 1 



X 

-1 



with total optical depth tjv along line-of-sight A: 
Qn 



TN = 



-7"0.66 



0.66 



(13) 



(14) 
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and Tiq{y) the partial optical depth through the CSE be- 
tween us and the point y. The dust is assumed to be in 
radiative equilibrium with the diluted and attenuated stel- 
lar radiation field. The stellar radiation is generated by 
two stars that have temperatures T^, and (prT-^, and lumi- 
nosities and (jjhLi,. We adopt Q oc for the overall 
frequency dependence of the extinction coefficient, with p 
between ~ 1 and 2. The temperature profile is then given 
by (see also Sopka et al. 1985): 



7(^)^^+^(exp(^)-l)"'x 

X exp (^2??^ (5"-2d^) dr) (15) 
with u) = l/{4:+p) and 
A! 



4(3+p)!Z(4 + p) 



(16) 



where the values of the Euler Zeta function Z{x) = 
^J^^ for p = 1 and p = 2 are Z{5) 1.037 and 
Z{6) ~ 1.0173, respectively, and 
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he 



The correction factor 7 is: 



^ Mi 



exp(?7/r^) 



exp (77/ (</.tT^))-1 



(17) 



(18) 



where 7 = 1 in the case of a single star and 7 = 2 in the 
case of two identical stars. 

Model profiles are calculated for the mass-loss histories 
examined above. The profiles are convolved with a gaus- 
sian of (T = 0.26", which corresponds to the diffraction 
limit of the ESO 3.6 m at 10.1 /zm. We adopt a stellar 
temperature = 7500 K and (j)T = 4000/7500, with a 
luminosity ratio of 0^ = 1/1.8 (Molstcr et al. 1999, van 
Winckel et al. 1999). The only free parameter left is the 
stellar radius ii* of the primary. We fitted by eye, aim- 
ing at equal emission in the centre for all a. The results 
for p = 1 are plotted in Fig. 11 (bottom panel) and sum- 
marised in Table 1. 

4.3. Model results 

From the output parameters R, Rq, i?* and C, the mass-loss 
rates at the inner and outer radii (Mq and M), the total 
mass McsE in the dusty CSE and the total luminosity 
L of the object arc derived, assuming a distance D = I 
kpc. The results for p = 1 are summarised in Table 1. 
The emission profiles for p = 2 are indistinguishable from 
those for p = 1, with the only difference being the stellar 
radius (and hence the stellar luminosity) which is almost 
six times larger for p = 1. The CSE radii are related to 
the time in the past when the material at those radii had 
been expelled, assuming a constant expansion velocity of 



''^oxp = 35 km s~^. The stellar radius of the hot star is ex- 
pressed both in arcseconds and in solar radii. Again, these 
conversions are made under the assumption of a distance 
D = I kpc. Radii (not angular), dynamical times, and 
mass-loss rates scale linearly with distance, whereas the 
CSE mass and the (total) luminosity scale quadratically 
with distance. 

The mass-loss history (or alpha) is not much con- 
strained by the optical depth and the 10 /um emission. The 
dynamical times do not depend strongly on the mass-loss 
history, and we find almost identical total CSE masses 
of M = 0.027 MgD^fkpcJ^. Our model does not repro- 
duce the 10 /im emission exactly. Apart from an observed 
point-source component in the centre with a flux density of 
80 X 10"^^ W m-Vm"^ (= 90 Jy) that is not reproduced 
by our model, the observed emission is more compact than 
our model produces. Lowering the effective temperature of 
the star(s) and/or including absorption of stellar light by 
a fixed column density before the stellar radiation reaches 
the spatially resolved dust shell both lead to a larger de- 
rived stellar radius and luminosity — and thus reduces 
the need for an additional point source in the centre — 
but the shape of the emission profiles remains virtually 
identical. 

5. Discussion 

5.1. Binary orbit 

The preliminary binary velocity curve that we present can 
only be used as an indication for the maximum masses of 
the stellar components of the binary in AFGL 4106. Their 
current masses are limited to ^15sin~^i Mq, which is 
consistent both with the masses and luminosities derived 
for a distance of 3.3 kpc by Molster et al. (1999) and with 
the high ratio of 60 /um flux density over CO brightness 
temperature (see Josselin et al. 1996) if the inclination 
angle i is close to 90° (edge-on). 

5.2. Structure of the ionized nebula 

Our coronographic imaging and spectroscopy provide new 
insight into the complex structure of the circumstellar en- 
velope of AFGL 4106. Emission from Ha and the forbid- 
den lines [N 11] and [S 11] originates mainly from a thin 
shell envelopping a roughly spherical thicker shell of dust. 
This dust shell extends from ^ 2.7 x 10'^^D[kpc] cm to 
~ 6.6 X 10^^£'[kpc] cm and has an optical depth around 
Ha close to unity. There is evidence for part of the neb- 
ular emission to also arise from within the dust shell: the 
electron density as traced by the ratio of the two compo- 
nents of the [S 11] doublet is higher at smaller projected 
distances from the star. A mixed dust/ion shell may result 
from a clumpy dust distribution, with some of the stellar 
UV radiation able to permeate the dust shell. Indeed the 
10 ^m image of AFGL 4106 (Molster et al. 1999) sug- 
gests a clumpy density or temperature distribution. An 
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Table 1. Modelling parameters of the spatial profiles of the optical depth at 0.66 /im and the 10 fim flux density. 
Five models are concerned, with p = 1 and different mass-loss rate histories parameterised by a. Distance dependent 
quantities are given assuming a distance D = 1 kpc. 



# 


a 


R["] 






*o[yr] 


R4"] 






c 


M[M0 yr-i] 


Mo[M0 yr-i] 


M[Mq] 


itot [L0] 


A 


-2 


7 


948 


1.9 


257 


3.4x10" 


-4 


73 


0.05 


1.0 X lO"-' 


1.4 X 10"" 


0.027 


1.5x10* 


B 


-1 


4.9 


664 


1.85 


251 


3.1x10- 


-4 


67 


0.28 


4.1 X 10"^ 


1.1 X 10-* 


0.026 


1.3x10"* 


C 





4.4 


596 


1.8 


244 


3.0x10- 


-4 


64 


0.60 


7.9 X 10"® 


7.9 X 10"® 


0.028 


1.2x10* 


D 


+1 


4 


542 


1.7 


230 


2.9x10- 


-4 


62 


1.05 


1.3 X 10"" 


5.3 X lO"--' 


0.028 


1.1x10* 


E 


+2 


3.7 


501 


1.6 


217 


2.8x10- 


-4 


60 


1.60 


1.8 X 10-" 


3.3 X 10"^ 


0.027 


1.0x10* 



extinction map such as Ha/H/3 would be useful to trace 
the dust, whereas a radio map obtained in the free-free 
continuum may be used to trace the distribution of ion- 
ized matter. Density-bounded ionized nebulae like the one 
around AFGL 4106 are rare and presumably short-lived, 
but not unique: similar nebulae are observed around Lu- 
minous Blue Variables (e.g AG Car) and Wolf-Rayet stars 
(Voors et al., in preparation). 

Unless a significant population of very small grains is 
postulated, dust temperatures in the detached dust shell 
around AFGL 4106 do not reach more than a few x 10^ K 
and hence thermal emission from dust grains cannot ex- 
plain the part of the near-IR emission observed by Garcfa- 
Lario et al. (1994) that is (projected) co-spatial with the 
dust shell. Thus it seems more likely that the near-IR 
emission is due to either free-free emission from the ion- 
ized component in the dust shell or starlight that is scat- 
tered by the large grains (a ^ l/im) found by Molster 
et al. (1999). The line intensity ratios of Ha, [N 11] and 
[S 11] unambiguously point at excitation conditions simi- 
lar to those foimd in Planetary Nebulae (see Garcfa-Lario 
et al. 1991). In the bow-shaped region where the line emis- 
sion reaches its maximum the excitation mechanism may 
have a shock component, resulting in enhanced emission 
by [N 11] over that by Ha. The shock may be due to the 
collision of the CSE of AFGL 4106 with the interstellar 
medium if the system is moving with a supersonic veloc- 
ity with respect to the interstellar medium. Ha emission 
from regions outside of the bow shape (up to about twice 
the outer radius of the dust shell) may be explained by an 
H II region being formed. 

5.3. Evidence for on-going mass loss 

AFGL 4106 is experiencing on-going mass loss. Evidence 
for this is seen in enhanced absorption in the short- 
wavelength wings of Si 11, Fe i,ii and N i lines. Fluorescence 
emission of [Ca 11] (see also Riera et al. 1995) originates 
from material that has been expelled not more than a few 
decades ago. Furthermore, the core emission in the 10 /im 
image must be of circumstellar origin (either dust or ion- 
ized gas). If dust is present close to the star then it must 
be optically thin at visual wavelengths. The reason for this 
is that our modelling of the 10 fim emission yields a lu- 



minosity that is in perfect agreement with the integrated 
bolomctric luminosity and reasonable estimations for in- 
terstellar extinction (Garcia-Lario et al. 1994; Molster et 
al. 1999) if the dust extinction coefficient index p = 1 (as 
in Volk & Kwok 1988). 

5.4. Mass-loss history 

The mass contained in the dust shell is a few 
xlO-^£)^[kpc]^ M0 and ~ 13 times less than the mass 
as derived by Molster et al. (1999). This is in part due to 
the choices by Molster et al. for the (smaller) inner and 
(larger) outer radii of the dust shell that together yield a 
duration of the dusty mass loss t — that is about three 
times as long as we estimate here. They assume a gas-to- 
dust ratio of 100 whereas here we adopt a value of 200, but 
the tinic^-averaged dust mass-loss rate that Molster et al. 
derive from their modelling of the spectral energy distri- 
bution is an order of a magnitude higher than we estimate 
here. This may be the result of differences in the optical 
constants of the dust species used in their modelling using 
large grains (o ~ 1/um) from those of the dust in Volk & 
Kwok (1988) that we use here. The mass-loss rate estimate 
from the CO emission, however, is in between the mass- 
loss rates that we estimate for the inner and outer radii of 
the dusty CSE, and thus seems to favour our more mod- 
erate mass-loss rates over the extremely high mass-loss 
rates estimated by Molster et al. (1999). The discrepancy 
might be resolved by storing most of the dust mass into 
large grains at large distances from the star, outside the 
CSE as traced by the 10 /xm emission, yielding the ob- 
served emission at wavelengths longward of ^ 50/im. 

The expansion velocity has increased over time from ~ 
30 km s"^ at the time when the material at the outer rim 
of the dust shell was expelled (tracer: [N 11]) through ^ 35 
km s"^ (tracer: CO) up to ~ 40 km s"^ at present (tracers: 
Si II, Fe I, II and N i). In combination with a steady mass 
loss this velocity evolution could yield the radial density 
dependence of pocr"^'^ as derived by Molster et al. (1999). 
This consideration suggests that of our models, model C 
is to be favoured: a constant mass-loss rate during the 
formation of the dust shell. 
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